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Abstract 

The  Rio  de  la  Plata,  one  of  the  largest  rivers  on  Earth,  discharges  into  the  ocean  waters  from  basin  that 
covers  a  large  area  of  South  America.  Its  plume  extends  along  northern  Argentina,  Uruguay,  and  southern 
Brazil  shelves  strongly  influencing  the  ecosystems.  In  spite  of  this,  little  is  known  about  the  mechanisms 
that  control  it.  Here  we  report  results  of  simulations  with  POM  carried  out  to  investigate  the  roles  of  wind 
and  river  discharge  in  Plata  plume  dynamics.  Different  outflows  were  explored,  including  an  average 
climatological  value  and  magnitudes  representative  of  La  Nina  and  El  Nino.  Forcing  the  model  with  river 
discharge  the  average  plume  speed  was  directly  related  to  the  outflow  intensity.  The  Plata  northward 
extension  varied  from  850  to  1550  km  and  for  average  discharge  a  band  of  low  salinity  waters  formed 
from  the  estuary  up  to  30°N  of  South  Brazilian  Shelf.  Upwelling  and  downwelling  winds  were  applied 
after  130  days.  The  distribution  of  low  salinity  waters  over  the  shelf  was  more  sensitive  to  the  wind 
direction  than  to  the  river  outflow  variability.  Downwelling  winds  were  very  capable  of  advecting  the  low 
salinity  signal  downshelf.  Upwelling  winds  were  efficient  in  eroding  the  plume,  which  was  basically 
detached  from  the  coast  by  Ekman  drift.  Abnormal  plume  intrusions  toward  low  latitudes  may  be  a  result 
of  the  original  plume  position  coupled  with  events  of  persistent  strong  downwelling  favorable  winds. 
Descriptors:  Eastern  South  America  Continental  Shelf,  South  Brazilian  Shelf,  Rio  de  la  Plata,  South  Brazil 
Bight,  Buoyancy  Plume,  Princeton  Ocean  Model. 

Resumo 

O  Rio  da  Prata,  um  dos  maiores  rios  da  Terra,  descarrega  no  oceano  aguas  de  uma  bacia  de  drenagem  que 
cobre  uma  ampla  area  da  America  do  Sul.  Sua  pluma  extende-se  ao  longo  do  norte  da  Argentina,  Uruguay 
e  sul  do  Brasil  influenciando  amplamente  os  ecossistemas  costeiros.  A  despeito  disso,  pouco  se  sabe  a 
respeito  dos  mecanismos  que  a  controlam.  Relatamos  aqui  simuIagBes  conduzidas  com  o  modelo  POM  na 
investigafao  do  papel  dos  ventos  e  da  descarga  fluvial  na  dinamica  da  pluma  do  Prata.  Descargas  com 
valores  mddios  climatologicos  e  magnitudes  representativas  de  El  Nino  e  La  Nina  foram  explorados. 
Somente  sob  descarga  fluvial,  a  velocidade  mddia  de  penetragao  da  pluma  ao  longo  da  costa  foi 
diretamente  relacionada  &  vazSo.  Sua  extensSo  variou  entre  850  e  1550  km  e  no  caso  da  vazSo  m6dia  uma 
banda  de  baixa  salinidade  formou-se  desde  o  estuario  at6  30°N  na  plataforma  continental  sul  brasileira.  O 
efeito  de  ventos  de  sudocste,  que  causam  subsidencia,  e  de  nordeste,  que  promovem  a  ressurggneia 
costeira,  foram  investigados  apos  130  dias  de  descarga  fluvial.  Os  resultados  sugerem  que  a  distribuigao  de 
baixa  salinidade  6  muito  mais  sensfvel  a  direfao  dos  ventos  do  que  &  descarga  fluvial.  Ventos  de  sudoeste 
sSo  capazes  de  advectar  o  sinal  de  baixa  salinidade  ao  longo  da  costa.  Ja  os  ventos  de  nordeste 
demonstram  ser  eficientes  na  erosao  da  pluma,  que  6  destacada  da  costa  por  deriva  de  Ekman.  Uma 
intrusSo  anormal  da  pluma  de  baixa  salinidade  em  direjSo  a  baixas  latitudes  pode  ser  o  resultado  da 
posipao  original  da  pluma  acoplada  a  eventos  persistentes  de  fortes  ventos  de  sudoeste.  Descritores: 
Plataforma  continental  leste  Sul  Americana,  Plataforma  continental  Sul  Brasileira,  Bacia  de  Santos,  Rio  da 
Prata,  pluma  costeira,  modelagem  hidrodinamica. 

Descriptors:  Eastern  South  America  Continental  Shelf,  South  Brazilian  Shelf,  Rio  de  la  Plata,  South 
Brazil  Bight,  Buoyancy  Plume,  Princeton  Ocean  Model. 
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Introduction 

The  intermittent  intrusion  of  relatively  fresh 
cold  waters  from  higher  latitude  into  the  South  Brazil 
Bight  (SBB),  between  28°S  and  22°S  off  the  Brazilian 
coast,  is  a  phenomenon  that  has  intrigued  the  scientific 
community.  This  phenomenon  has  been  known  for  a 
long  time  but  was  first  documented  in  the  scientific 
literature  in  the  middle  90’s  (Campos  et  al,  1995, 
1996a,b).  Combining  in  situ  measurements  and 
satellite  imagery,  these  studies  demonstrated  the 
presence  of  cold  (T~18°C)  low  salinity  (S~32)  waters 
extending  from  south  of  the  Cape  of  Santa  Marta 
(28°S)  to  the  middle  of  the  bight,  at  approximately 
23°30’S  (Fig.  1).  The  absence  of  major  river  systems 
in  the  region  and  the  thermohaline  characteristics  of 
the  coastal  waters  excluded  the  possibility  of  influence 
of  local  river  runoff  or  upwelling  (Campos  et  al, 
1996b). 

Based  on  satellite  SST  images  and 
trajectories  of  surface  drifters  launched  in  February 
and  late  April  of  1993,  the  unexpected  water  mass  was 
associated  with  advection  from  the  vicinity  of  the 
Brazil-Malvinas  Confluence  (BMC)  in  the  Argentinian 
continental  shelf  (Stevenson  et  al.,  1998;  Campos  et 
al,  1996b).  This  water  mass  was  suggested  as  the 
agent  responsible  for  the  failure  of  important  Brazilian 
fisheries  (Matsuura,  1996;  Sunye  &  Servain,  1998) 
and  for  the  intrusion  of  organisms  from  different 
biogeographic  zones.  The  detection  of  cold  water 
foraminifera  at  unexpected  low  latitudes  (Stevenson  et 
al.,  1998)  and  the  anecdotal  reports  of  the  occasional 


appearance  of  penguins  are  some  indicators.  It’s 
recurrence  was  studied  by  Campos  et  al.  (1999)  and 
Lentini  et  al.  (2001)  through  the  analysis  of  a  13-year 
time  series  of  AVHRR  images.  Their  findings  revealed 
the  cold  intrusions  as  seasonal  phenomenon,  which 
presents  an  interannual  variability  possibly  related  to 
El  Nino/Southem  Oscillation  (ENSO).  At  the  same 
time,  Piola  et  al.  (2000)  based  on  the  analysis  of 
historical  hydrographic  data  described  the  presence  of 
a  defined  seasonal  sign  associated  with  the  Rio  de  la 
Plata  surface  plume.  As  one  of  the  largest  rivers  on 
Earth,  the  Plata  River  (Plata  hereafter)  discharges  into 
the  ocean  waters  from  a  drainage  basin  that  covers 
nearly  20%  of  South  America.  The  influence  of  its 
plume,  as  suggested  by  the  historical  data,  extends  into 
the  continental  shelf  off  northern  Argentina,  Uruguay 
and  southern  Brazil.  Sometimes  the  plume  reaches 
latitudes  as  low  as  27°S  during  the  winter,  while  in  the 
summer  it  is  confined  to  32°S  (Fig.  2)  (Piola  et  al., 
1999,  2000,  2005). 

The  Plata  discharge  presents  interannual 
variability  related  to  El  Nino  scale  rainfall  variations 
over  South  America  (Ropelewski  &  Halpert,  1987; 
Grimm  el  al.,  1998;  Mechoso  &  Iribarren,  1992; 
Depetris  et  al.,  1996).  This  correlation  was  suggested 
as  a  reason  for  unusual  plume  intrusions  along  the 
shelf,  which  under  high  outflow  events  would  explain 
the  appearance  of  cold  and  fresh  waters  to  SBB 
(Campos  et  al.,  1999  and  Piola  et  al.,  2000).  Surface 
salinity  observations  from  Guerrero  et  al.  (1997a)  in 
late  April  and  May  of  1983  demonstrated  how  the  30 
isohaline  extended  to  about  150  km  offshore  and  to 


|sallnilY/.  winter 


SST  /  winter  1993 


longitude 


Fig.  1.  a)  Sea  surface  temperature  (SST)  derived  from  AVHRR/NOAA  for  July  1993.  The  dark  colors  indicate  cold  waters 
coming  from  the  southwest  into  the  area,  b)  Surface  salinity  field  derived  from  COROAS-WOCE  hydrographic  1993 
winter  cruise  (Modified  from  Campos  et  al,  1999). 
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200  km  along  the  coast  of  Uruguay  beyond  its  average 
position  during  a  high  outflow  event.  The  1993 
intrusion  however,  remains  unclear.  The  change  in  the 
extent  was  of  the  order  of  800  to  1000  km,  or  nearly  4 
times  larger  than  observed  by  Guerrero  et  at.  (1997a). 
Nevertheless,  the  anomalous  northward  penetration 
occurred  about  one  year  after  a  large  river  outflow 
event  associated  with  the  1992-1993  El  Nifio. 

Despite  its  considerable  dimensions  and 
strong  influence  on  nearshore  ecosystems  (Ciotti  et  al., 
1995;  Mianzan  et  al.,  2001),  little  is  known  about  the 
Plata  plume  variability  or  mechanisms  that  control  it 
A  few  numerical  modeling  studies  investigated  its 
seasonal  variability  (Simionato  et  al.,  2001),  most  of 
the  investigations  focused  on  the  estuarine  domain 
(O’Connor,  1991;  Framinan  &  Brown,  1996;  Guerrero 


et  al.,  1997a, b;  Simionato  et  al.,  2004).  The 
phenomenon  of  coastal  buoyant  plumes,  on  the  other 
hand,  is  well  explored  in  the  scientific  literature. 
Works  related  to  the  Amazon  River  plume  (Lentz, 
1995;  Limebumer  et  al.,  1995;  Geyer  et  al.,  1996),  the 
Delaware  and  Chesapeake  systems  (MOnchow  and 
Garvine,  1993;  Rennie  et  al,  1999;  Sanders  and 
Garvine,  2001;  Houghton  et  al.,  2004),  Columbia 
(Hickey  et  al.,  1998)  and  the  Rhine  river  plumes 
(Simpson  and  Souza,  1995)  are  some  examples.  These 
studies  reveal  the  strong  dependence  that  buoyant 
plumes  demonstrate  not  only  under  river,  but  under 
wind  forcing  as  well  (Chao  &  Boicourt,  1986;  Chao, 
1987,  1988a,b;  Kourafalou  et  al.,  1996a,b;  Garvine, 
1999). 


XL  (Km) 


Fig.  2:  a)  Orthogonal  rotated  bathymetric  domain,  extending  from  south  of  Rio  de  la  Plata  estuary  to 
Cabo  Frio,  b)  Zoom  at  estuarine  zone,  c)  Estuarine  cross-section  and  sigma  layers. 
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Within  this  context,  the  objective  of  the 
present  study  is  to  extend  the  understanding  on  the 
Plata  plume  variability  and  investigate  the  relative 
roles  of  river  discharge  and  wind  forcing  on  the 
dynamics  of  the  Plata  river  plume  over  the 
Southwestern  Atlantic  shelf.  Specifically,  we  strove  to 
explain  plume  intrusions  to  abnormal  low  latitudes. 
This  research  was  based  on  simple  numerical 
simulations  of  the  Plata  plume  using  the  Princeton 
Ocean  Model  (POM). 

Methodology 

Buoyancy  driven  coastal  plumes  have  been 
widely  investigated  by  many  researchers.  These 
features  are  nonlinear,  with  Coriolis  force  and  friction 
playing  important  roles  (Csanady,  1984;  Munchow  & 
Garvine,  1993;  Kourafalou  et  al.,  1996  b,a;  Garvine, 
1999).  Here  we  adopted  the  Princeton  Ocean  Model 
(POM),  a  three-dimensional,  sigma-coordinate, 
primitive  equation  model  (Mellor,  1996;  Blumberg  & 
Mellor,  1987). 


The  Numerical  Model 

The  model  dependent  variables  are  the 
zonal,  meridional  and  upward  components  of  velocity 
(u,  v,  w,  respectively),  the  free  surface  elevation  q, 
temperature  T,  salinity  S,  and  two  properties  of  the 
turbulence  field  (turbulent  macro  scale  and  kinetic 
energy).  The  density  p  is  computed  from  the  equation 
of  state.  The  system  of  equations  include  the  nonlinear 
advective  terms  and  the  variables  are  computed 
integrating  the  continuity  equation,  the  momentum 
equation,  and  two  equations  for  advection-diffusion  of 
heat  and  salt.  For  the  horizontal  grid,  an  “Arakawa  C” 
differencing  scheme  is  used.  The  horizontal  time 
differencing  scheme  is  explicit,  whereas  the  vertical 
differencing  is  implicit,  which  eliminates  time 
constraints  for  the  vertical  coordinate  and  permits  the 
use  of  fine  vertical  resolution  in  surface  and  bottom 
layers.  An  imbedded  second  moment  turbulent  closure 
sub-model  (Mellor  &  Yamada,  1982)  provides  the 
vertical  mixing  coefficients,  while  the  horizontal 
diffusion  and  viscosity  coefficients  are  parameterized 
by  the  Smagorinsky  formula  (Blumberg  &  Mellor, 
1987). 


set  to  dx=dy=10  km,  which  resulted  in  a  105x270 
points  grid  (Fig.  2a).  With  this  spatial  resolution,  the 
Nyquist  wavelength  is  of  20km,  a  bit  larger  than  the 
local  Rossby  Radius  of  deformation.  However, 
considering  both  the  alongshelf  and  cross-shelf 
dimensions  of  the  plume,  the  10km  horizontal 
resolution  is  reasonable.  Depths  greater  than  1000  m 
were  set  to  1000  m  and  the  vertical  structure  was 
represented  by  12  sigma  surfaces  (Fig.  2c).  In  the 
numerical  simulations,  the  domain  was  initialized  with 
a  homogeneous  structure,  with  T=25°C  and  S=35  psu. 
The  Plata  River  outflow  was  implemented  as  an  open 
boundary  condition  at  the  estuary  head  on  the  western 
side  of  the  domain.  Since  the  domain  was  rotated  40° 
degrees,  the  estuary  head  fit  well  to  the  meridional 
domain  axis  (Fig.  2b).  Thus,  the  horizontal  barotropic 
velocities  across  the  river  mouth  could  be  prescribed 
as: 

U=~,  y=0,  (1) 

A 


where  A  is  the  area  of  the  vertical  section  of  the 
estuary  head,  and  Q  is  the  river  discharge  [m3,s-1].  The 
free  surface  elevation  is  computed  by  the  model, 
applying  a  no-gradient  condition  (qB-i  =qB,  where  B 
is  the  boundary  index).  Hydrographic  data  show  that 
the  water  column  at  the  estuary  head  is  quasi 
homogeneous  in  temperature  and  salinity  (Guerrero  et 
al.,  1997a,b).  Temperature  and  salinity  were  specified 
using  a  temporal  relaxation  based  on  historical 
conditions  of  the  Rio  de  la  Plata  winter  temperature 
(T=14°C)  and  salinity  (S=5  psu). 

Oceanic  Boundary  Conditions 


In  addition  to  the  river  runoff  specification 
as  an  open  boundary  condition  at  the  estuary  head, 
three  oceanic  open  boundary  conditions  were 
implemented.  At  the  cross-shore  boundaries  (north  and 
south)  the  elevation  was  treated  with  a  no  gradient 
condition,  while  barotropic  velocity  components  and 
the  horizontal  internal  velocity  u  and  v  were  calculated 
by  the  Sommerfeld  radiation  condition: 


dr  ~dr 
— +C— 
dt  dy 


=  0, 


(2) 


Domain  Setting,  Initial  Conditions  and  River  Forcing 
Implementation 

The  model  was  implemented  in  a  physical 
domain  that  encompasses  the  southeastern  South 
American  continental  shelf  from  40°S  to  22°S.  The 
model  domain  is  a  regular  orthogonal  grid  rotated  40° 
from  the  north.  The  bathymetric  data  used  in  its 
construction  were  based  on  Brazilian  Navy  nautical 
charts,  Argentinean  Hydrographic  Naval  Service  and 
the  ETOP05  database.  The  horizontal  resolution  was 


where  <{>  denotes  the  velocity  component,  C  is  the 
phase  velocity,  and  y  is  the  axis  normal  to  the  open 
boundary.  The  evaluation  of  C  usually  follows  the 
shallow  water  wave  velocity  expression,  but  in  this 
case,  it  was  evaluated  by  a  simple  diagnostic  version 
proposed  by  Orlanski  (1976): 
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Here,  the  superscript  denotes  the  time  step,  and  the 
subscript  the  position  related  to  the  boundary.  Once  C 
is  evaluated,  the  velocity  component  can  be  calculated 
at  the  next  time  step  n+1  by: 


r 


n+1  _ 
B  ~ 


rr,(1-r)+2r;„ 

1+r 


with 


1  if 
C  if 
0  if 


Cal 

0<C<1 

C<0 


(4) 


At  the  eastern  oceanic  boundary,  a  similar 
implementation  to  the  cross-shore  boundaries  was 
performed,  but  with  the  addition  of  a  10-point 
hyperbolic  sponge  layer  for  the  velocity  components. 
In  that  layer,  the  velocity  components  were 
progressively  damped  from  their  original  value  (4>b-io) 
inside  of  the  domain,  to  78%  of  their  value  near  the 
boundary  (0.78<(>b-i)-  An  Orlanski  scheme  was  applied 
at  the  eastern  domain  edge  (<|>B).  For  the  temperature 
and  salinity  boundary  conditions,  an  advection 
solution  was  used  in  order  to  allow  the  plume 
movement  through  the  cross-shore  boundaries: 


dS  dS 

- +V - 

dt  dy 


=  0, 


(5) 


where  v  is  the  velocity  component  normal  to  the 
boundary,  and  S/Sy  is  the  derivative  in  the  normal 
direction.  This  equation  was  implemented  through  a 
scheme  advanced  in  time  and  space,  as  described  by 
Palma  &  Matano  (2000): 

^=S’B+ri(S'B-S'B^)+r2(S'^-S'B)  (6) 

where 

r=0.5(d/,/dy)(^+|YnB|) 

r2=0.5(dr/dy)(v"s-|v£|). 

In  these  equations,  S  is  the  salinity,  Ati  is  the 
internal  mode  time  step.  Ay  is  the  meridional  space 
resolution,  vB  is  the  normal  velocity  component  of  the 
boundary  and  Sc  is  the  prescribed  external  value. 
Temperature  is  treated  in  an  identical  fashion.  On  the 
eastern  boundary,  a  simple  no-gradient  condition  was 
applied  for  temperature  and  salinity,  since  it  was  far 
enough  to  ensure  no  influence  on  the  plume. 


Numerical  Experiments 


In  the  real  ocean,  estuarine  plume  dynamics 
are  dictated  by  several  factors.  The  river  discharge 
interacts  with  the  surrounding  water,  generating 
thermohaline  gradients  that  drive  the  density 
circulation;  the  wind  works  by  establishing  sub- 
inertial  coastal  currents;  the  tides  enhance  vertical 
mixing;  and  nearby  oceanic  currents  impose  horizontal 


shear.  Here  we  focused  on  only  two  forcing  factors: 
the  river  outflow  and  the  wind  forcing.  To  understand 
the  role  of  these  mechanisms  and  to  facilitate  the 
interpretation  the  model  runs  were  divided  into  two 
sets  of  experiments.  First,  we  conducted  simulations 
with  different  river  outflows,  and  them  combined  river 
outflow  with  different  scenarios  of  wind  forcing. 

Experiments  Forced  Only  with  River  Outflow 

In  this  set  of  experiments,  we  dealt  with  the 
buoyancy-driven  force  generated  by  the  Plata 
discharge,  which  mixes  with  the  surrounding  oceanic 
waters.  Our  experiments  were  based  on  two  major 
questions.  (1)  What  would  be  the  extent  of  the 
alongshelf  penetration  of  the  plume  for  a  steady 
coastal  buoyancy  source?  (2)  How  would  the 
alongshelf  penetration  vary  with  increased  or 
decreased  river  discharges? 

To  address  these  questions,  we  first  ran  three 
experiments  with  the  following  values  of  river 
discharge:  Q=15xl03mV  (low);  Q=25xl03m3s-1 
(medium);  and  Q=50xl03m3s_l  (high).  According  to 
Mechoso  and  Iribarren  (1992)  and  Depetris  et  al. 
(1996)  these  values  are  typical  of  a  La  Nifia  situation 
(dry  season/lower  discharge),  to  normal  periods 
(average  discharge),  and  to  El  Nino  Events  (higher 
discharge).  Hereafter,  we  refer  to  these  experiments  as 
R15,  R25  and  R50.  The  river  outflow  was  set  to  zero 
at  initialization  and  allowed  to  grow  linearly  to  the 
target  value  within  one  inertial  period.  The 
experiments  were  run  for  lengths  of  time  enough  for 
allowing  the  progression  of  the  plume  along  the  entire 
extent  of  the  spatial  domain.  The  results  were  analyzed 
by  observing  the  spatial  evolution  of  the  34.9 
isohaline,  as  depicted  in  Figure  3.  Within  the  first  20 
days  of  integration,  the  plume  remained  confined  to 
the  estuarine  domain,  except  for  the  higher  outflow 
case  (R50),  when  the  34.9  isohaline  extended  beyond 
Punta  del  Este  during  that  period.  In  all  cases,  after 
reaching  the  shelf  region,  the  plume  turned 
anticlockwise,  following  the  direction  of  propagation 
of  coastally  trapped  waves.  From  that  point  on,  the 
main  difference  among  the  three  experiments  was  the 
northward  extension  of  the  plume.  In  the  R15 
experiment,  after  80  days  of  simulation  with  the  low 
river  discharge  the  plume  extremity  reached  latitude 
34°S,  an  extent  of  about  623  km.  In  the  R25  and  R50 
experiments,  the  distance  from  the  plume  extremity  to 
the  estuary  head  was  764  km  and  1099  km, 
respectively.  On  day  140,  the  R15  plume  had  reached 
the  latitude  of  Rio  Grande  (32°S),  while  the  R25 
plume  extended  to  the  northern  limit  of  Patos  Lagoon 
(30°S).  In  R50,  the  plume  reached  27°S,  past  the  Cape 
of  Santa  Marta  Grande.  The  alongshelf  plume 
penetration,  as  determined  by  the  distance  of  the  34.9 
isohaline  from  the  estuary  head,  is  represented  in 
Figure  4a.  Within  the  first  200  km,  inside  the  estuary, 
the  propagation  velocities  of  the  plumes  in  the  three 
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Fig.  3:  Surface  salinity  field  evolution  for  R15,  R25  and  R50  river  only  experiments.  The 

contours  depict  the  34.9  isohaline  position,  ranging  from  day  20  to  140  with  30-day  increments. 

experiments  are  more  or  less  similar,  suggesting  an  increase  in  river  discharge.  This  can  be  seen  in  the  first 
estuarine  channel  flux.  When  the  plumes  reach  the  derivative  for  the  velocity  curve,  depicted  by  the  solid 

distance  of  400  km  from  the  estuary  head,  their  line  in  Figure  4b.  An  upshelf  plume  penetration 

velocities  become  notably  dissimilar.  This  leads  to  an  (southwestward)  was  also  noticeable  in  these 

increased  distance  between  their  frontal  positions.  To  experiments,  achieving  significant  distances  (more 

verify  the  behavior  of  the  plume  as  a  function  of  river  than  250  km)  along  Argentinean  coast  (Fig.  3).  The 

discharge  in  a  more  continuous  fashion,  we  ran  distance  of  this  intrusion  was  also  proportional  to  the 

another  set  of  experiments,  changing  the  river  out-flow  river  discharge  magnitude.  None  of  these  plumes  fully 

intensity  by  smaller  increments.  The  average  speed  of  achieves  an  alongshore  stationary  position  of  the 
propagation  of  the  plume  head,  as  a  function  of  river  plume  head  during  the  simulations.  However,  with 

outflow,  can  be  followed  in  the  starred  graph  in  Figure  respect  to  offshore  expansion,  the  plumes  seem  to 

4b;  the  speed  of  propagation  increases  as  discharge  achieve  a  limit.  A  considerable  change  in  plume  cross¬ 
increases.  In  our  main  experiments,  the  estimated  shore  dimensions  was  noticed  with  increasing  river 

values  were  5.88  km. day"1  (0.07  m.s~‘)  for  R1 5,  7.84  discharge,  but  all  cases  eventually  displayed  a 

km.day  1  (0.09  m.s"1)  for  R25,  and  11.35  km. day-1  tendency  toward  a  cessation  of  the  offshore 

(0.13  m.s  f)  for  R50.  We  note  that  the  rate  of  increase  development.  The  differences  in  offshore  extension 

of  the  velocity  decreases  with  increased  discharge,  can  be  verified  by  plotting  plume  cross-shore  width  as 
which  suggests  that  for  higher  outflow  regimes  there  is  a  function  of  time  along  two  cross-shore  transects,  one 
a  weaker  response  (in  terms  of  plume  velocity)  to  an  at  the  estuary  mouth,  the  other  off  the  latitude  of 
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Lagoa  Mirim  (transects  J=34  and  J=80,  in  Figure  5c). 
The  result  is  shown  in  Figure  4c.  On  the  northern 
transect  (J  =  80),  the  34.9  isohaline  distance  from  the 
coast  reaches  60  km  for  R15,  70  km  for  R25,  and  110 
km  to  R50.  Along  the  southern  transect  (J  =  34),  the 
cross-shore  expansion  stabilizes  at  approximately 
360  km  for  R15,  420  km  for  R25,  and  500  km  for 
R50.  Figure  5a  shows  the  surface  salinity  distribution 
for  the  experiment  with  average  river  discharge 
(R25).  That  figure  reveals  the  freshening  of  a  region 
extending  from  the  northern  Argentinian  coast  to 
southern  Brazil,  in  the  form  of  a  low  salinity  band. 
The  main  difference  between  this  spatial  distribution 
and  the  higher  discharge  case  (not  shown)  was  a 
further  alongshore  extension,  with  a  more 
pronounced  freshening  of  waters  along  the  Brazilian 
inner  shelf,  and  a  wider  offshore  dislocation  of  the 
estuarine  front.  Figure  5  also  shows  the  horizontal 
salinity  distribution  near  the  bottom  (5b)  and  cross 
sections  along  three  cross-shelf  transects  (5c).  In 
spite  of  the  wider  cross-shore  expansion  of  the  plume 
at  the  ocean  surface,  the  salinity  field  near  the 
bottom  was  roughly  confined  by  the  40-meter 
isobath.  A  similar  situation  was  verified  in 
experiments  R15  and  R50  (not  shown  here).  The 
vertical  distribution  was  stratified  with  positive 
slope  of  the  salinity  isopleths,  as  seen  in  the  three 
transects  J=34,  J=80  and  J=122  (Fig.  6c).  The 
estuarine  salinity  configuration  is  evident  along  J  =  34 
where  the  34.9  isohaline  extends  450  km  from  the 
estuary  head.  To  the  north  of  the  estuary  mouth,  the 
isohaline  configuration  is  consistent  with  an 
alongshore  flux  (northeastward)  induced  by 
geostrophic  adjustment. 

Near  the  estuary  head,  the  flow  is  nearly 
straight,  following  the  main  channel.  The  flow 
regime  changes  when  low-salinity  water  reaches  the 
sea  and  starts  to  accumulate,  forming  a  large  bulge. 
This  region  displays  a  complicated  meandering 
pattern,  showing  vortices  with  diameters  of 
approximately  35  km  and  velocities  of  about  0.125 
m.s-1,  as  indicated  by  transect  J  =  34  in  Figure  6.  The 
transition  between  this  meandering  field  and  the 
coastal  body  is  smooth,  within  200  km  north  of  the 
estuary  axis.  In  this  zone,  a  surface  coastal  jet 
develops.  At  the  core  of  the  J  =  80  transect, 
northeastward  alongshore  velocities  up  to  0.2  m.s”1  are 
found.  This  velocity  decreases  to  the  north,  decaying 
to  less  than  0.1  m.s-1  across  transect  J=122. 
Underneath  the  surface  jet,  a  weak  counter  current  is 
observed  near  the  bottom  on  transect  J  =  80.  It  is 
mainly  composed  of  slightly  diluted  near-bottom 
waters,  and  -0.025  m.s”1  speed  southwestward.  In 
simulations  R15  and  R50  the  circulation  followed 
almost  the  same  general  pattern,  differing  only  in 
terms  of  magnitude,  being  0.13  m.s”1  and  0.25  m.s”1 
respectively,  for  surface  velocities  (across  J  =  80),  and 
-0.023  m.s”1  and  -0.035  m.s”1,  respectively  for  bottom 
velocities  (J  =  80). 


Fig.  4.  a)  Alongshelf  plume  penetration  from  the  estuary  head 
as  a  function  of  time  for  all  numerical  experiments,  b) 
Average  velocity  of  plume  propagation  as  function  of  river 
discharge  only  (denoted  by  the  stars).  The  thin  dashed  line 
represents  a  regression  to  the  data  in  the  form  of  a  power 
curve  (std=0.0022;  R2=0.9954),  while  the  solid  line 
represents  its  first  derivative,  (c)  Cross-shore  plume  width  for 
river  only  experiments  and  simulations  under  upwelling- 
favorable  winds  (the  cross-sections  are  indicated  on  Figure 
6). 

Experiments  Forced  With  Wind  and  River  Discharge 

In  this  second  group  of  simulations, 
initialization  files  from  the  river  forcing  experiments 
were  used  to  start  the  runs.  For  this  new  set,  we 
selected  all  three  (R15,  R25  and  R50)  discharge 
experiments.  These  new  simulations  started  on  day 
130.  The  river  discharge  was  kept  constant  since  the 
initialization,  while  the  wind  forcing  was  ramped  up 
within  one  inertial  period.  After  initialization,  the 
winds  were  kept  constant  in  time  and  throughout  the 
domain.  The  plume  response  to  different  alongshore 
wind  directions  was  tested.  Since  the  grid  was  rotated 
40°  with  respect  to  the  north,  we  assumed  NE  and  SW 
winds  to  be  coincident  with  the  domain’s  meridional 
axis.  The  main  wind  experiments  were  forced  with  8 
m.s”1  wind  speed,  for  both  northeasterly  and 
southwesterly  cases.  This  corresponds  to 
|ry|=10”4xm2.s”2,  where  ry  is  the  wind  stress  divided  by 
the  water  density.  For  convenience,  we  will  refer  to 
these  experiments  as  R15NE,  R25NE  and  R50NE, 
respectively,  for  low,  medium  and  high  discharge 
northeasterly  wind  simulations,  and  R15SW,  R25SW 
and  R50SW,  respectively,  for  low,  medium  and  high 
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Fig.  5:  Surface  (a)  and  bottom  (b)  salinity  fields  for  R25  river  only  experiment  after 
140  days  of  simulation,  c)  salinity  transects  for  the  same  experiment.  The  40  m  isobath 
is  depicted  by  a  dashed  line  on  “b”.  Salinity  is  depicted  by  the  contours  10,  20,  30,  32, 
33,  34  and  34.9  PSU,  as  well  by  a  shade  from  fresh  (black)  to  salt  (white)  waters. 


Fig.  6.  Meridional  velocity  field  for  R25  river  only 
experiment  at  day  140.  Contour  intervals  are  between  - 
0.2  and  0.2  m.s'1  with  0.025  increment.  Negative 
velocities  are  depicted  by  dashed  lines. 
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discharge  southwesterly  wind  simulations.  In  the 
southwesterly  wind  case,  two  other  wind  magnitudes 
were  used  to  force  the  model,  referred  to  as  weaker  (6 
m.s”'~Ty=-5x  lCf5m2.s-2)  and  stronger  (15 
m.s”1~Ty=-3x  10”4m2.s”2)  southwesterly  winds 
experiments.  The  tx  component  was  set  to  zero  in  all 
simulations. 

Southwesterly  Wind  Experiments 

In  this  set  of  experiments,  we  applied  the  so- 
called  downwelling  winds  over  low,  medium  and  high 
discharge  plumes  for  50  days.  We  will  first  describe 
the  evolution  of  the  34.9  surface  isohaline  along  the 
coast  as  a  function  of  time  for  the  main  wind 
experiment  (Fig.  7).  We  will  then  focus  on  the 
differences  imposed  by  wind  magnitude.  As  an  initial 
condition,  the  plume  in  the  low  discharge  experiment 
(R15SW)  extends  850  km  from  the  estuary  to  Rio 
Grande  (32°S);  in  the  average  discharge  experiment 
(R25SW)  the  plume  head  is  near  Tramandal  (30°S), 
approximately  1100  km  from  the  mouth;  and  in  the 


high  discharge  experiment  (R50SW)  the  head  is 
located  about  450  km  farther  upcoast,  near  27°S. 

Ten  days  after  southwesterly  winds  were 
applied,  the  plume  in  all  three  experiments  extended 
approximately  240  km  northeastward  of  their  initial 
positions,  reaching  Tramandai  for  R15SW,  Santa 
Marta  Cape  (28.3°S)  for  R25SW  and  25.2°S  for 
R50SW  (near  Paranagua).  A  significant  modification 
in  the  offshore  expansion  of  the  plume’s  bulge  was 
noted,  with  a  contraction  of  the  low  salinity  region 
towards  the  estuary.  After  20  days,  the  high  discharge 
plume  shifted  about  275  km  farther  up-coast  and 
reaches  the  middle  portion  of  the  South  Brazil  Bight 
at  23.8°S.  In  R25SW,  the  plume  head  was  located 
at  26.2°S,  450  km  downcoast,  and  on  R15SW  it  was 
found  at  the  latitude  of  Santa  Marta  Cape.  After  30 
days  of  southwesterly  wind  forcing,  the  plume  in  the 
R25SW  experiment  entered  SBB,  east  of  ParanaguA, 
while 'on  R50SW  experiment  the  plume  occupied 
almost  all  SBB  domain.  The  plume  head  reached 
the  middle  of  the  bight  in  R25SW  and  R15SW, 
after  40  and  50  days  of  simulation  respectively.  These 


Fig.  7.  Surface  salinity  field  evolution  for  R15SW,  R25SW  and  R50SW 
downwelling-favorable  wind  experiments  (8  m.s"1  magnitude).  The  contours  depict 
the  34.9  isohaline  position  with  5-day  increments. 
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differences  in  plume  extension  with  and  without  winds 
are  summarized  in  Figure  8.  After  day  130,  the 
application  of  downwelling  winds  led  to  a  substantial 
modification  of  the  plume’s  propagation  speed  for  all 
river  discharge  conditions,  with  an  average  velocity  of 
27. 16  km. day”1  (0.31  m.s-1,  about  2  to  4.5  times  faster 
than  the  river  only  cases).  In  an  attempt  to  evaluate  the 
influence  of  wind  strength  in  alongshore  velocity  of 
the  plume,  two  additional  downwelling  experiments 
were  performed  for  all  river  outflows,  with  stronger 
and  weaker  southwesterly  winds  (15  m.s-1  and  6 
m.s’1).  The  plume  propagation  speed  was  21.93 
km. day’1  (0.25  m.s’')  and  42.22  km.day’1  (0.49 
m.s’1),  respectively,  for  the  weaker  and  stronger  wind 
experiments.  A  summary  of  the  results  in  all 
southwesterly  cases  is  shown  in  Table  1.  The  plume 
propagation  speed  was  independent  of  the  river 
outflow.  In  Figure  4a  we  observe  the  along  coast 
evolution  in  these  two  additional  wind  experiments  for 
an  average  river  outflow  magnitude.  The  dashed  line 
represents  stronger  (denoted  by  the  numeral  I)  and 
weaker  (denoted  by  II)  wind  cases.  Since  the 
propagation  speed  was  independent  of  river  outflow 
magnitude,  similar  paths  for  the  other  experiments  can 
be  found  by  drawing  parallel  lines  from  the  end  of  the 
simulation  forced  only  with  river  outflow.  A  surface 
salinity  field,  when  the  plume  reaches  the  S.  Sebastiao 
latitude,  is  presented  in  Figure  8  for  the  main  wind 
experiments  (8  m.s’1  southerly  winds),  where  R25SW 
displays  its  40th  day  of  simulation  (top  panel)  while 
R50SW  is  in  its  25lh  day  (bottom  panel).  Both  plumes 
extended  more  than  2000  km  from  the  Plata  estuary. 
Although  representing  different  situations,  these 
experiments  were  very  similar  in  a  general  way.  In  its 
central  part,  between  the  Uruguay  coast  and 
Florian6polis,  the  low  salinity  band  was  composed 
mainly  by  20  to  30  PSU  waters  that  displayed  cross¬ 
shore  salinity  gradients.  Salinities  below  20  PSU 
remained  mostly  confined  to  the  estuary,  whereas 
higher  alongshore  salinity  gradients  were  observed 
near  the  head  of  the  plume.  Small  low  salinity  pools 
were  frequently  detached  at  the  estuary  and  near  the 
Uruguay  coast  in  these  simulations.  Other 
downwelling  wind  experiments  produced  only  slightly 
different  surface-salinity  configurations.  In  Figure  9 
we  observe  R25SW  salinity  transects  after  40  days  of 
southwesterly  winds.  The  inclined  configuration  of 
isohalines  observed  in  the  outflow  only  experiment 
was  replaced  by  a  cross-shore  salinity  gradient,  with  a 
lack  of  vertical  salinity  stratification  (see  vertical 
isohalines  in  transects  J=80,J=122  and  J=202).  The 
plume  widened  beyond  the  60  m  isobath.  At  the 
estuary  domain  (J=34)  there  was  a  retraction  of  the 
plume  body  and  an  enhancement  of  the  salinity 
gradients.  This  configuration  was  different  from  the 
outflow  only  simulations,  where  the  front  occupied  the 
whole  estuary.  In  transects  J=34  and  J=80,  we  observe 
surface  buoyant  pools  detached  from  the  plume’s  main 
body,  with  a  length  scale  of  30  km. 


Fig.  8.  Surface  salinity  fields  for  downwelling- 
favorable  wind  experiments  with  average  discharge 
outflow  after  40  days  of  wind  forcing  (top  panel)  and 
extreme  river  discharge  after  25  days  of  wind  (bottom 
panel).  Salinity  is  depicted  by  the  contours  of  10,  20, 
30,  32,  33,  34  and  34.9  PSU  as  well  by  the  shade 
ranging  from  fresh  (black)  to  salt  (white)  waters. 


Fig.  9.  Salinity  transects  for  R25SW  after  40  days  of 
wind  forcing  with  8  m.s'1  strength. 
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This  rapid  extension  of  the  low-salinity 
signal  along  the  coast  is  a  result  of  advection  promoted 
by  a  barotropic  flux  over  the  shelf.  The  downwelling 
wind  induced  a  surface  Ekman  transport  towards  the 
coast,  accumulating  water  in  the  nearshore  region. 
This  was  enhanced  near  the  estuarine  zone,  where 
higher  surface  elevations  were  observed  against  the 
Uruguayan  coast,  achieving  magnitudes  of  0.25  m, 
0.40  m  and  0.8  m,  respectively,  for  the  weak,  medium 
and  strong  wind  cases.  Along  the  southern  Brazilian 
coast,  the  elevation  varied  between  0.15  and  0.5  m, 
while  along  the  SBB  the  largest  elevation  did  not 
exceed  0.45  m  for  the  stronger  wind  case.  The 
elevation  pattern  followed  closely  the  depths  over  this 
region  for  all  cases,  revealing  the  confinement  of  the 
flux  on  the  shelf. 

This  piling  up  of  water  creates  a  positive 
cross-shore  pressure  gradient  which  is  governed  by 
Coriolis  adjustment  producing  a  strong  and  well 
defined  northeastward  current.  In  Figure  10  we 
observe  the  alongshore  velocity  component  for  four 
transects  in  the  domain  after  20  days  of  SW  wind 


forcing.  A  definite  barotropic  signal  is  found  on  the 
continental  shelf,  nearly  bounded  by  the  250-300  m 
isobaths.  This  current  started  on  the  southern 
boundary,  traveled  through  the  north  Argentinian, 
Uruguayan,  and  Brazilian  shelves.  A  greater 
intensification  occurred  on  transects  J  =  80  and  J  = 
122  for  south  Brazilian  shelf.  Currents  achieved 
velocities  of  0.8  m.s”1  and  0.7  m.s”1,  respectively,  in 
the  core,  but  mostly  reached  the  average  value  of  0.1- 
0.3  m.s-1  on  middle  and  outer  shelves.  Wider  currents 
in  the  northern  and  southern  shelf  transects  achieved 
0.45  m.s”1  and  0.65  m.s”1,  respectively,  in  their  core. 
Transect  J=202  had  the  smallest  alongshore  velocity 
gradient  across  the  shelf.  In  other  wind  experiments, 
currents  achieved  maximum  values  of  1.0  m.s”1  and 
0.55  m.s”1,  respectively,  for  the  weak  and  strong  wind 
cases.  In  the  estuarine  zone,  the  current  followed  the 
bathymetry,  making  a  wide  turn  and  entering  the 
region  immediately  outside  of  the  Plata  estuary.  This 
current  shear  was  responsible  for  the  detachment  and 
advection  of  small  low-salinity  pools  (Fig.  9,  J=34). 
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Northeasterly  Wind  Experiments 

Northeasterly  (upwelling  favorable)  winds 
were  applied  over  developed  plumes  as  in  the  previous 
experiments,  but  for  6  m.s-1  and  8  m.s-1  wind  values. 
The  model  was  run  for  20  days.  The  surface-salinity 
fields  for  days  5  and  10  are  shown  in  Figure  11  for 
average  and  high  outflow  cases  under  8  m.s-1  NE 
winds.  The  plume  responded  to  wind  forcing  within  a 
few  days  with  a  significant  offshore  spread.  Within  5 
days,  a  considerable  offshore  dislocation  was 
produced.  For  both  the  river  discharge  cases  the 
plumes  extended  45  km  offshore.  After  10  days  of 
wind  forcing,  the  coastal  plume  achieved  almost  twice 
its  original  cross-shore  width.  The  plume  separation 
from  the  coast  increased  at  the  plume  nose  and  on  the 
southern  Uruguayan  shelf.  For  the  high  discharge  case, 
the  separation  also  occurred  on  southern  Brazilian 
coast  east  of  Rio  Grande.  On  day  15  (not  shown)  the 


plume  detachment  continued,  and  at  the  end  of  the 
simulation  its  body  was  completely  separated  from  the 
coast.  The  plume  offshore  drift  may  be  seen  by  the 
evolution  of  the  low-salinity  signal  in  the  estuary 
mouth  and  Lagoa  Mirim  salinity  transects  on  R25NE 
(Fig.  4c)  for  both  wind  magnitudes.  The  plume  grew 
offshore  in  an  approximately  linear  fashion,  doubling 
its  width  on  the  northern  transect  (J=80)  after  10  days 
of  wind  forcing.  The  average  growth  was  similar  for 
other  river  outflows  on  the  estuarine  transect:  8 
km. day-1  for  6  m.s-1  winds  (dashed  lines  with 
numerals  III  and  IV)  and  near  15  km. day-1  for  8  m.s-1 
winds  (continuous  lines).  The  velocities  were  the  same 
for  the  high  outflow  case  on  J=80,  although  for  the 
average  and  lower  outflow  cases  they  diminished  after 
some  time  (i.e.,  once  the  plume  nose  started  to  get 
closer  to  the  transect  position).  In  terms  of  along  coast 
retraction,  the  modification  is  clear  in  all  NEwind 
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Fig.  11.  Surface  salinity  field  after  5  and  10  days  under  upwelling-favorable  wind  forcing  with  average 
(a,b)  and  extreme  outflow  (c,d).  The  dashed  line  represents  the  initial  34.9  isohaline  position. 
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Fig.  12.  Salinity  transects  forR25NE  after  15  days  of  wind 
forcing. 


current,  this  retraction  could  reach  approximately  -18 
km. day-1  for  8  m.s-1  winds  and  -12  km. day-1  for  6 
m.s”1  winds,  on  average  and  high  outflow  cases  during 
the  first  10  days.  But  once  these  plume  bodies 
departed  from  the  coast  and  from  the  shelf  current, 
their  southward  rate  of  propagation  was  diminished. 
For  example,  we  can  follow  the  alongshore  path 
evolution  measured  by  the  model  for  the  average  river 
discharge  case  under  8  m.s”1  \vinds  (Fig.  4a,  depicted 
as  R25NE).  The  difference  between  river  outflow-only 
and  the  northeastward  experiment  regimes  was  noticed 
by  the  curve  break  on  this  graph.  About  10  days  affer 
northeasterly  winds  were  applied  the  main  plume  body 
was  detached  100  km  from  the  coast  (Figs  lib  and 
lid),  and  only  a  narrow  low-salinity  strip  remained 
attached  to  the  coast  (Fig.  12).  The  main  plume  body 
was  10  to  15  m  thick,  while  the  coastal  strip  extended 
throughout  the  water  column  between  the  coast  and 
the  10  m  isobath.  In  the  estuary  domain  the 
northeasterly  winds  dislocated  the  plume  offshore  in 
the  first  5  days.  After  that,  a  southwestward  current 
developed  and  reversed  the  plume.  As  exemplified  by 
10  and  20  isohalines  (Fig.  11),  it  changed  its  natural 
up-coast  propagation  for  a  down-coast  dislocation. 
This  shelf  current  occurred  as  a  result  of  the 
geostrophic  balance  between  the  negative  pressure 
gradient  (developed  by  the  offshore  surface  Ekman 
transport)  and  the  Coriolis  force.  After  20  days,  for 
both  weaker  and  stronger  winds,  the  near-coast  water 
level  decreased  by  -0.20  and  -0.55  m  at  the  estuary  and 


by  -0.20  and  -0.30  m  on  the  south  Brazilian  shelf.  The 
currents  were  about  -0.35  m.s-1  for  stronger  winds  (8 
m.s”1),  and  -0.20  m.s”1  for  weaker  winds  (6  m.s”1). 
The  barotropic  current  follows  the  bathymetry 
southwestward,  and  flows  through  the  southern 
boundary.  Northeasterly  winds  led  to  significant  and 
rapid  extension  of  low  salinity  water  along  the  coast 
south  of  the  estuary. 

Discussion 

Fed  by  the  second  largest  river  system  of 
South  America  with  an  average  river  runoff  of 
23xl03m3.s”'  the  Plata  shelf  plume’s  dimension 
surpass  the  internal  Rossby  radius  and  therefore  is 
subject  to  earth  rotation  effects.  The  funnel  shape  of 
the  estuary  provides  a  gradual  change  between  an 
estuarine  channel  and  the  inner  shelf,  with  an 
increasing  Kelvin  number  (K  =  b/ri,  where  b=inlet 
breadth;  ri=  internal  Rossby  radius  of  deformation  on 
the  shelf)  (Garvine,  1999).  We  observed  distinctive 
circulation  types  along  the  estuarine  axis  that  agree 
with  the  regime  differentiation  found  in  the  numerical 
studies  of  Chao  &  Boicourt  (1986).  These  features 
were  also  qualitatively  consistent  to  observations  of 
Miinchow  &  Garvine  (1993)  for  the  Delaware  plume, 
despite  the  differences  in  the  spatial  scale  compared  to 
the  Plata.  Inside  the  estuarine  channel  (first  200  km), 
we  found  velocities  consistently  higher  than  those  over 
the  shelf.  After  leaving  this  “channel  regime”,  the 
water  flowed  out  of  the  estuary  and  traveled  up-coast, 
but  a  considerable  amount  of  energy  seems  to  be 
accumulated  near  the  estuary  mouth,  forming  a  large 
pool  of  mixed  water.  A  “near  flow  field”  characterized 
the  anticyclonic  surface  flow  and  the  formation  of  the 
estuarine  bulge.  The  “far  field”  was  characterized  by 
the  plume  alongshelf  intrusion. 

Near/far  field  separation  can  be 
characterized  by  a  strong  cyclonic  surface  flow  and 
downwelling  (Chao  &  Boicourt,  1986;  Oey  &  Mellor, 
1993),  and  with  a  modification  from  a  nearly  inertial 
regime  in  the  turning  region,  to  a  nearly  geostrophic 
regime  in  the  coastal  current  itself  (Garvine,  1987; 
Miinchow  &  Garvine,  1993).  Although  no  sharp 
distinction  was  observed  in  our  simulations,  the 
differentiation  between  near  and  far  fields  was 
apparent  in  river  only  experiments.  In  spite  of  this 
regime  differentiation,  the  bulge  surface  area  was 
related  to  outflow  magnitude,  and  a  meandering 
pattern  was  a  common  feature  in  this  location. 

The  presence  of  upshelf  intrusions  or 
massive  bulges  such  as  those  observed  in  our 
simulations  are  usually  absent  in  field  observations 
(Garvine,  2001).  One  exception  is  the  upshelf 
intrusion  observed  of  the  Changjiang  (Yangtze)  river 
plume,  as  described  by  (Beardsley  et  al.,  1985,  in 
Garvine,  2001).  Fong  &  Geyer  (2002)  attributed  these 
differences  to  the  effect  of  ambient  shelf  currents, 
which  would  be  essential  in  the  determination  of  the 
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plume  structure  by  models.  Garvine  (2001),  on  the 
other  hand  demonstrated  that  simulated  plumes  are 
highly  sensitive  to  model  configuration.  Significant 
changes  on  bulge  or  in  upshelf  penetration  may  result 
from  different  schemes  of  outflow  implementation. 
Despite  of  these  findings,  the  formation  of  a  bulge  and 
the  upshelf  intrusion  for  the  Plata  simulations  were 
consistent  to  historical  hydrographic  maps  of  Piola  et 
al  (2000),  and  the  salinity  distributions  of  Guerrero  et 
al.  (1997a). 

In  terms  of  salinity  distribution  for  the 
estuary,  our  simulations  resembled  quite  well  Guerrero 
el  al.  (1997a)  and  Piola  et  al.  (2000)  observations, 
where  10-20  PSU  isohalines  are  located  near  Punta 
Rasa  and  Punta  del  Este.  The  vertical  stratification  and 
plume  bottom  confinement  could  also  be  reasonably 
reproduced,  matching  Guerrero  et  al.  (1997a)  mild 
wind  synoptic  situation.  Our  experiments  also 
demonstrated  a  small  but  noticeable  offshore  shift  of 
the  estuarine  front  with  increase  of  river  outflow. 
Following  an  alongshore  path  toward  northeast,  with  a 
coastal  length  varying  from  850  km  to  1550  km 
(respectively  for  low  and  high  outflow  only  cases),  the 
Plata  exerted  strong  influence  over  the  Uruguayan  and 
South  Brazilian  shelves.  The  salinity  fields  in  our  river 
outflow  only  experiments  corroborate  Miranda  (1972) 
analysis  of  coastal  water  mass  formation,  where  a  low 
salinity  band  near  the  South  Brazil  Shelf  (SBS)  was 
associated  with  Plata  waters. 

These  results  also  agreed  with  the  known 
seasonal  salinity  variability  of  the  Plata  plume  (Piola 
el  al.,  2000).  In  summer,  under  the  average  discharge 
of  Q  =  20*  lOW1,  the  plume  mean  location  stays 
near  Rio  Grande,  while  during  the  winter 
(30x1 0!mV)  its  influence  extends  to  27°S  on  the 
SBS.  These  positions  depart  to  some  extent  from  those 
in  our  outflow  cases,  but  most  of  the  differences  may 
be  attributed  to  wind  seasonality. 

Simpson  &  Souza  (1995)  studying  the  action 
of  tides  upon  a  buoyancy-driven  current  demonstrated 
its  effect  on  vertical  mixing  through  elevated  bottom 
stress  and  tidal  straining  of  the  stratified  flow 
generated  by  the  plume.  Garvine  (1999)  investigated 
the  effect  of  tidal  amplitudes  on  the  alongshelf 
penetration  of  a  coastal  buoyant  plume  in  an  idealized 
domain.  A  linear,  inverse  relationship  suggested  that 
regions  dominated  by  extreme  tidal  amplitude  (e.g. 
Bay  of  Fundy)  are  expected  to  have  greatly  reduced 
alongshelf  penetrations.  On  the  other  hand,  in 
continental  shelves  such  as  SBS  where  tides  have 
small  amplitude,  plume  intrusions  are  expected  to  be 
largely  unaffected  by  them.  In  addition,  Simionato  et 
al.  (2001)  numerical  simulations  demonstrated  that 
tides  may  play  a  role  in  the  Plata  plume  mixing 
process,  but  mainly  by  smoothing  the  horizontal 
salinity  structures. 

The  northeastward  geostrophic  flow  along 
SBS  was  consistent  with  the  modeling  performed  by 
Pereira  (1989)  and  Ghisolfi  &  McKee  (2003)  under 


thermohaline  forcing.  Our  results  also  reinforce  the 
findings  of  a  buoyancy-driven  current  described  by 
Zavialov  et  al.  (1998)  by  means  of  a  inverse  model  for 
circulation  and  by  Soares  &  Mdller  (2001)  who 
computed  northward  residual  flow  from  a  current 
meter  off  Rio  Grande  shelf.  The  weaker  counter- 
current  underneath  the  northeastward  surface  flow  was 
also  described  in  Ghisolfi  &  McKee  (2003).  These 
resultant  “two  layer  circulation”  flows  are  a  common 
feature  in  coastal  plumes  and  may  be  recognized  as  the 
first  mode  of  a  baroclinic  Kelvin  wave  (Chao,  1987; 
Garvine,  1999). 

The  plume’s  alongshelf  average  speed  of 
propagation  was  linked  to  the  outflow  magnitude,  but 
faster  growth  was  observed  for  the  lower  river  outflow 
magnitude.  Apparently,  the  water  accumulation  at  the 
bulge  and  plume  cross-width  growth  compensate  for 
the  lack  of  increase  in  alongshore  growth  rate.  We 
verified  that  even  with  high  outflow,  the  river 
simulations  could  not  explain  penetrations  beyond 
27°S.  This  model  behavior  motivated  the  exploration 
of  wind  forced  simulations  since  winds  play  an 
important  role  in  SBS  currents.  To  study  the  plume 
response  to  steady  winds,  we  carried  out  simulations 
that  included  alongshore  winds  under  different  river 
discharges.  Downwelling  winds  induced  an  onshore 
plume  contraction  due  to  flux  towards  the  coast.  The 
water  piling  up  near  the  coast  leads  to  the 
establishment  of  an  offshore  pressure  gradient  and  a 
northeastward  current.  This  barotropic  flow  was 
responsible  for  the  alongshore  plume  growth.  Our 
results  show  that  the  mean  intrusion  velocities  under 
wind  forcing  were  independent  of  the  river  outflow 
magnitude.  The  effect  of  downwelling-favorable 
winds  has  been  investigated  by  Chao  (1987)  using  an 
idealized  domain,  Kourafalou  et  al.  (1996b)  for  South 
Atlantic  Bight  (U.S.A),  Hickey  et  al.  (1998)  for 
Columbia  &  Rennie  et  al.  (1999)  for  Chesapeake 
plume.  Their  results  displayed  similar  effects  coherent 
to  our  observations:  a  reduction  of  the  plume’s  surface 
area  and  a  downcoast  acceleration  in  a  jet-like 
configuration. 

The  resulting  velocities  over  the  shelf  were 
consistent  with  the  findings  of  Zavialov  et  al.  (2002) 
on  the  SBS,  who  measured  0.6-0.8  m.sH  currents 
under  strong  alongshore  winds.  In  some  cases,  the 
observed  currents  achieved  1.2  m.s-1  at  15  m  depth. 
Under  persistent  positive  alongshore  currents  Zavialov 
et  al.  (2002)  recorded  a  reduction  of  surface  salinity 
from  35  to  30  PSU.  These  values  were  maintained  for 
more  than  30  days  in  late  autumn.  Upwelling- 
favorable  winds  exerted  very  different  effects. 
Regardless  of  the  wind  and  outflow  intensities,  surface 
Ekman  drift  acted  to  separate  the  plume’s  body  from 
the  coast.  Previous  studies  have  demonstrated  this 
response  by  means  of  field  observations  (Miinchow  & 
Garvine,  1993;  Sanders  &  Garvine,  2001;  Houghton  et 
al,  2004)  and  numerical  modeling  (Chao,  1987; 
Kourafalou  etal.,  1996b;  Fong  &  Gey er,  2001).  The 
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offshore  plume  dislocation  were  comparable  to 
findings  of  Houghton  et  al.  (2004)  for  Delaware 
plume,  The  detachment  from  the  coast  agrees  well 
with  their  observations  and  with  the  conceptual  model 
proposed  by  Fong  &  Geyer  (2001).  The  persistence  of 
winds  with  a  strong  upwelling  component  are  likely  to 
be  an  effective  mechanism  for  dissipating  the  Plata 
plume  and  transferring  shelf  waters  to  the  ocean. 
Kourafalou  et  al.  (1996a)  have  shown  how  upwelling 
winds  work  on  the  exchange  of  materials  from  the 
continental  shelf  and  Piola  (2002)  suggested  this 
transfer  of  light  water  pools  to  the  Brazil  Malvinas 
Confluence  as  an  effective  and  irreversible  process  of 
water  mass  exchange.  Near  the  coast,  simultaneously 
to  the  plume  detachment  occurred  the  upwelling  of 
subsurface  waters.  This  process  seemed  to  be 
enhanced  off  the  Uruguayan  coast,  where  the  effect  of 
winds  upon  the  shelf  currents  was  also  capable  of 
reversing  the  river  plume  direction.  This  result  was 
consistent  with  Mianzan  et  al.,  (2001)  observations  in 
the  summer  of  2000.  In  a  period  of  anomalously  low 
Plata  outflow  (12xl03mV )  and  strong  northeasterly 
winds,  high  salinity  waters  (>  34)  and  Southern 
Brazilian  marine  biota  (plankton  and  fishes)  were 
observed  in  the  outer  Plata  estuary. 

In  general,  our  results  suggest  that  the 
distribution  of  low-salinity  waters  over  the  shelf  is 
much  more  sensitive  to  wind  direction  than  to  the  river 
outflow  variability.  Alongshore  winds  seem  to 
determine  the  main  role  on  the  Plata  plume  alongshore 
growth  or  destruction.  The  effect  of  cross-shore  wind 
components,  as  investigated  in  Pimenta  (2001)  have 
shown  smaller  importance  on  the  plume  control.  Once 
intrusion  velocities  did  not  depend  on  the  river 
discharge  magnitude  but  mostly  on  the  downwelling 
wind  intensity,  the  possible  intrusion  of  cold  and  fresh 
waters  up  to  the  SBB  seems  to  be  linked  to  three  main 
factors:  the  initial  alongshelf  plume  position  (that 
might  be  determined  by  the  river  outflow),  the 
intensity,  and  the  persistence  of  the  southwesterly 
wind  events.  Our  simulations  have  show  that  the  time 
scale  necessary  to  advect  the  Plata  waters  to  SBB  by 
steady  wind  forcing  and  average  river  discharge, 
varied  between  25  and  49  days,  depending  on  the  wind 
strength  (Table  1).  When  constant,  strong 
downwelling  winds  blow  over  a  plume,  with  spatial 
extensions  similar  to  the  mean  winter  field  described 
by  Piola  et  al.  (2000),  this  time  scale  can  be  even 
reduced  to  15  days. 


An  important  part  of  the  shelf  currents 
variance  along  the  Southeastern  America  continental 
shelf  is  known  to  be  attributed  to  the  winds  (Soares  & 
Moller,  2001;  Zavialov  et  al.,  2002;  Pimenta  et  a!., 
2004).  Cold  fronts  are  responsible  for  strong 
southwesterly  winds  and  for  the  generation  of  sea  level 
oscillations  that  travels  equatorward  along  the  SBS 
and  SBB  in  the  form  of  continental  shelf  waves  (Stech 
8c  Lorenzzetti,  1992;  Castro  &  Lee,  1995).  Unusual 
strong  events  might  occur.  Melo  et  al.  (2003)  describe 
an  cyclone  generated  at  Uruguay  as  responsible  for 
raising  sea  level  about  1  m  and  for  generation  of  fairly 
strong  northward  currents  (1.6  m.s-1)  at  the  Santa 
Catarina  shelf.  Indirect  drift  measurements  registered 
around  27.7°S  have  also  revealed  persistent  northward 
currents  for  10  to  almost  20  days  in  early  May  and 
June  2001  (Pimenta  et  al.,  2004).  Assireu  et  al.  (2003) 
suggest  that  the  drifter  that  accompanied  the  cold 
water  intrusion  of  1993  winter  (described  in  Campos 
et  al,  1996a)  was  also  driven  by  continental  shelf 
waves. 

Rather  than  being  the  result  of  continuous 
wind  forcing,  plume  intrusions  might  result  from  the 
passage  of  a  sequence  of  cold  fronts.  Rennie  et  al. 
(1999)  have  demonstrated  the  effect  of  winds  in  the 
creation  of  a  pulsed  buoyancy  current  downstream  of 
Chesapeake  Bay  &  Hickey  et  al.  (1998)  have  shown 
sensitive  responses  of  the  Columbia  plume  to  changes 
in  wind  speed  and  direction  in  the  time  scale  of  2-10 
days. 

Since  the  Plata  outflow  was  normal  during 
the  winter  of  1993,  the  observations  of  anomalously 
cold,  low  salinity  waters  are  probably  attributed  to 
southwesterly  winds.  Studies  demonstrated  up  to  a 
50%  increase  in  the  formation  of  cyclones  over 
Uruguay  and  San  Matias  Gulf  (Argentina)  during 
ENSO  events  (Satyamurty  et  al,  1990;  Gan  &  Rao, 
1991).  This  conclusion  is  in  agreement  with  recent 
analysis  of  historical  hydrographic  data  and  winds, 
which  suggest  that  the  winter  1993  southwesterlies 
were  among  the  strongest  since  1949  (Piola  et  al, 
2005).  Extreme  events  may  arise  from  certain 
combinations  of  alongshore  winds  and  outflow. 
Extreme  northeastward  penetrations  of  the  Plata  plume 
will  be  associated  with  large  outflow  and  strong, 
sustained  southeasterly  winds,  while  extreme  plume 
retractions  will  occur  under  low  outflow  and 
northeasterly  winds. 


Table  1.  Necessary  time  under  downwelling  winds  for  the  plume  to  reach 
S.  SebastiSo  (t,,)  and  alongshore  average  speed  of  propagation  (Vp)  under 
different  wind  and  river  magnitudes. 


Q  (xloW.s'1)  W  (days) 

_ i,=-5xl0'imJs'>  T,=-10J,mi.5'!  v^xlO’W.s'* 

15  >60  48.5  30.7 

25  49.7  39.9  25.5 

50  28.2  23.1  15.1 

V„  (km.day'1)-^  -21.93 _ -27,16 _ -42.22 
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Conclusions 

Numerical  simulations,  using  the  Princeton 
Ocean  Model,  revealed  the  effects  of  river  discharge 
and  wind  forcing  on  the  behavior  of  the  Rio  de  la  Plata 
plume  along  the  continental  shelf  of  southeastern 
South  America,  The  results  show  that  increased  river 
outflows  produce  deeper  northward  penetration  of  the 
plume  along  the  coast.  However,  in  the  absence  of 
favorable  winds,  the  time  scales  associated  with  the 
alongshore  penetration  of  the  plume  are  unrealistically 
long  to  explain  intrusions  in  the  SBB.  On  the  other 
hand,  the  experiments  with  the  inclusion  of  wind 
forcing  reveal  that  the  behavior  of  the  plume  is 
significantly  more  dependent  on  the  wind  pattern  than 
on  the  magnitude  of  the  river  discharge.  When  the 
winds  blow  from  the  southwest,  in  a  downwelling 
favorably  configuration,  the  penetration  of  the  plume 
along  the  coast  happens  in  a  very  coherent  way.  From 
our  results  we  infer  that  a  combination  of  the 
persistence  and  strength  of  southwesterly  winds  might 
lead  to  very  intense  events,  enhancing  the  plume 
penetration  along  the  coast,  which  could  explain  the 
episodes  of  anomalous  penetration  of  the  plume  into 
the  South  Brazil  Bight. 
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